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Abstract 

First, we analyze trajectories inside the Feigenbaum attractor and 
obtain the atypical weak sensitivity to initial conditions and loss of 
information associated to their dynamics. We identify the Mori sin- 
gularities in its Lyapunov spectrum with the appearance of a special 
value for the entropic index q of the Tsallis statistics. Secondly, the dy- 
namics of iterates at the noise-perturbed transition to chaos is shown 
to exhibit the characteristic elements of the glass transition, e.g. two- 
step relaxation, aging, subdiffusion and arrest. The properties of the 
bifurcation gap induced by the noise are seen to be comparable to 
those of a supercooled liquid above a glass transition temperature. 
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1 Introduction 



Evidence for the incidence of nonextensive dynamical properties at critical 
attractors in low dimensional nonlinear maps has accumulated and advanced 
over the last few years; specially with regards to the onset of chaos in logistic 
maps - the Feigenbaum attractor [1, 7], and at the accompanying pitchfork 
and tangent bifurcations [8, 9]. The more general chaotic attractors with 
positive Lyapunov coefficients have full-grown phase-space ergodic and mix- 
ing properties, and their dynamics is compatible with the Boltzmann-Gibbs 
(BG) statistics. As a difference, critical attractors have vanishing Lyapunov 
coefficients, exhibit memory-retentive nonmixing properties, and are there- 
fore to be considered outside BG statistics. 

Naturally, some basic questions about the understanding of the dynamics 
at critical attractors are of current interest. We mention the following: Why 
do the anomalous sensitivity to initial conditions and its matching Pesin 
identity obey the expressions suggested by the nonextensive formalism? How 
does the value of the entropic index q arise? Or is there a preferred set of 
q values? Does this index, or indexes, point to some specific observable 
properties at the critical attractor? 

From a broader point of view it is of interest to know if the anomalous 
dynamics found for critical attractors bears some correlation with the dy- 
namical behavior at extremal or transitional states in systems with many 
degrees of freedom. Two specific suggestions have been recently advanced, 
in one case the dynamics at the onset of chaos has been demonstrated to be 
closely analogous to the glassy dynamics observed in supercooled molecular 
liquids [10], and in the second case the dynamics at the tangent bifurcation 
has been shown to be related to that at thermal critical states [11]. 

With regard to the above comments here we briefiy recount the following 
developments: 

i) The finding [7] that the dynamics at the onset of chaos is made up of an 
infinite family of Mori's g-phase transitions [12, 13], each associated to orbits 
that have common starting and finishing positions located at specific regions 
of the attractor. Every one of these transitions is related to a discontinuity 
in the a function of 'diameter ratios' [14], and this in turn implies a g - 
exponential and a spectrum of g-Lyapimov coefficients equal to the Tsallis 
rate of entropy production for each set of attractor regions. The transitions 
come in pairs with conjugate indexes q and Q = 2 — q, as these correspond 
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to switching starting and finishing orbital positions. The amphtude of the 
discontinuities in a diminishes rapidly and consideration only of its dominant 
one, associated to the most crowded and sparse regions of the attractor, 
provides a very reasonable description of the dynamics, consistent with that 
found in earlier studies [1, 4]. 

ii) The realization [10] that the dynamics at the noisc-pcrturbcd edge of 
chaos in logistic maps is analogous to that observed in supercooled liquids 
close to vitrification. Four major features of glassy dynamics in structural 
glass formers, two-step relaxation, aging, a relationship between relaxation 
time and configurational entropy, and evolution from diffusive to subdiffusive 
behavior and finally arrest, are shown to be displayed by the properties of 
orbits with vanishing Lyapunov coefficient. The previously known properties 
in control-parameter space of the noise-induced bifurcation gap [14, 15] play 
a central role in determining the characteristics of dynamical relaxation at 
the chaos threshold. 

2 Mori's g-phase transitions at onset of chaos 

The dynamics at the chaos threshold [i = [icoi the 2;-logistic map 

j^{x)^l-ii\x\\ z>\-l<x<\ (1) 

has been analyzed recently [3, 7]. The orbit with initial condition = 

(or cquivalently, Xq = 1) consists of positions ordered as intertwined power 
laws that asymptotically reproduce the entire period-doubling cascade that 
occurs for < This orbit is the last of the so-called 'superstable' periodic 
orbits at JL^ < \Xc-,n— 1, 2, ... [14], a superstable orbit of period 2°°. There, 
the ordinary Lyapunov coefficient Ai vanishes and instead a spectrum of q- 
Lyapunov coefficients A^^^ develops. This spectrum originally studied in Refs. 
[13] when ^ = 2, has been shown [4, 7] to be associated to a sensitivity to 
initial conditions (defined as ^tixo) = limAa;o^o(Aa;t/Aa;o) where Aa;o is 
the initial separation of two orbits and Axt that at time t) that obeys the 
gr-exponential form 

Uxo) = expJA,(a;o)t] = [I - {q - l)\{xo) r'^'"' (2) 

suggested by the Tsallis statistics. Notably, the appearance of a specific value 
for the q index (and actually also that for its conjugate value Q = 2 — q) 
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works out [7] to be due to the occurrence of Mori's 'g-phase transitions' [12] 
between 'local attractor structures' at /Xc- 
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Figure 1: Absolute values of positions in logarithmic scales of iterations r for 

a trajectory at 

He with initial condition xq = 0. The numbers correspond to iteration times. 



As shown in Fig. 1, the absolute values for the positions x^- of the tra- 
jectory with Xt=o = at time-shifted t — t + 1 have a structure consist- 
ing of subsequences with a common power-law decay of the form r~^/^~^ 
with q = 1 — ln2/(2; — l)lnQ;(^) [3], where a{z) is the Feigenbaum uni- 
versal constant that measures the period-doubling amplification of iterate 
positions. That is, the attractor can be decomposed into position subse- 
quences generated by the time subsequences r = {2k + 1)2", each obtained 
by proceeding through n = 0, 1, 2, ... for a fixed value of A; = 0, 1, 2, .... See 
Fig. 1. The k = subsequence can be written as Xt = exp2_q(— A^^^t) with 
A(o) = {z-l)\na{z)/\n2. 

q-Lyapunov coefficients. The sensitivity Cti^o) can be obtained [7] from 
^t(m) ~ |(7„(m - l)/(7n(m)|", i = 2'^-!, n large, where (7n(m) = dn+i,m/dn,m 
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and where dn,m are the diameters that measure adjacent position distances 
that form the period-doubhng cascade sequence [14]. Above, the choices 
Axq — dn^rn ^nd Axt — dn^m+u t — 2"' — 1, hscvB been made for the initial 
and the final separation of the trajectories, respectively. In the large n limit 
(T„(m) develops discontinuities at each rational m/2"+^ [14], and according to 
our expression for ^i(m) the sensitivity is determined by these discontinuities. 
For each discontinuity of cr„(m) the sensitivity can be written in the forms 
— expg[Agt] and — ^W2-q[^2-qt], Xq > and A2-5 < [7]. This result 
reflects the multi-region nature of the multifractal attractor and the memory 
retention of these regions in the dynamics. The pair of g-exponentials corre- 
spond to a departing position in one region and arrival at a different region 
and vice versa, the trajectories expand in one sense and contract in the other. 
The largest discontinuity of (T„(m) at m = is associated to trajectories that 
start and finish at the most crowded (x ~ 1) and the most sparse {x ~ 0) 
regions of the attractor. In this case one obtains 

the positive branch of the Lyapunov spectrum, when the trajectories start 
at X ~ 1 and finish at x ~ 0. By inverting the situation one obtains 

g)^_ 2(.-l)lna(.) ^^0,1,2,..., (4) 
^ (2A; + l)ln2 ' ' ' ' ' w 

the negative branch of the Lyapunov spectrum. Notice that exp2_^(t/) = 
1/ expg(— y). So, when considering these two dominant families of orbits all 
the g-Lyapunov coefficients appear associated to only two specific values of 
the TsaUis index, q and Q — 2 — q. 

Mori's q-phase transitions. As a function of the running variable —00 < 
q < cxD the g-Lyapunov coefficients become a function A(q) with two steps 
located at q = g = 1 — In 2/(2; — 1) In a{z) and q = Q = 2 — g. In this manner 
contact can be established with the formalism developed by Mori and cowork- 
ers [12] and the g-phase transition obtained in Refs. [13]. The step function 
for A(q) can be integrated to obtain the spectrum 0(q) (A(q) = d(f)/dX{q)) 
and its Legendre transform ip{X) (= — (1 — q)A), the dynamic counterparts 
of the Renyi dimensions D{q) and the spectrum f{a) that characterize the 



5 



geometry of the attr actor. The result for ip{X) is 



(1 - Q)X, Ag) < A < 0, 
(1 -g)A, < A < Xf\ 



(5) 



As with ordinary thermal 1st order phase transitions, a " g-phase" transition 
is indicated by a section of linear slope m = 1 — g in the spectrum (free energy) 
'0(A), a discontinuity at q in the Lyapunov function (order parameter) A(q), 
and a divergence at q in the variance (susceptibility) ^(q)- For the onset of 
chaos at fic{z = 2) a g-phase transition was numerically determined [12, 13]. 
According to ip{X) above we obtain a conjugate pair of g-phase transitions 
that correspond to trajectories linking two regions of the attractor, the most 
crowded and most sparse. See Fig. 2. Details appear in Ref. [7]. 
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Figure 2: g-phase transitions with index values g = 0.2445 and Q = 2 — q = 
1.7555 obtained ior z — 2 from the main discontinuity in cr„(m). See text for 
details. 



Generalized Pesin identity. Ensembles of trajectories with starting points 
close to the attractor point Xq expand in such a way that a uniform dis- 
tribution of initial conditions remains uniform for all later times t. As a 
consequence of this we established [4, 7] the identity of the rate of entropy 
production K^''^ with A^^^. The g-generalized rate of entropy production Kg 



6 



is defined via Kqt = Sq{t) — Sq{0), t large, where 

(6) 

is the TsaUis entropy, Pi is the trajectories' distribution, and where In^ y = 
— — q) is the inverse of expg{y). See Figs. 2 and 3 in Ref. [4]. 

3 Glassy dynamics at noise-perturbed onset 
of chaos 

We describe now the effect of additive noise in the dynamics at the onset of 
chaos. The logistic map z — 2 reads now 

Xt+i = fnixt) = 1 - //Xj + Xtcr, -l<Xt<l, < // < 2, (7) 

where Xt is Gaussian-distributed with average {xtXf) = ^t.t', and a is the 
noise intensity. For a > the noise fluctuations wipe the fine features of the 
periodic attractors as these widen into bands similar to those in the chaotic 
attractors. nevertheless there remains a well-defined transition to chaos at 
/ic(o') where the Lyapunov exponent Ai changes sign. The period doubling of 
bands ends at a finite maximum period 2^^'^^ as — > /Xc(c) and then decreases 
at the other side of the transition. This effect displays scaling features and is 
referred to as the bifurcation gap [14, 15]. When a > the trajectories visit 
sequentially a set of 2" disjoint bands or segments leading to a cycle, but 
the behavior inside each band is fully chaotic. These trajectories represent 
ergodic states as the accessible positions have a fractal dimension equal to 
the dimension of phase space. When a — the trajectories correspond to 
a nonergodic state, since as t — > oo the positions form only a Cantor set of 
fractal dimension df = 0.5338.... Thus the removal of the noise a ^ leads 
to an ergodic to nonergodic transition in the map. 

As shown in Ref. [10] when /Xc(c > 0) there is a 'crossover' or 'relaxation' 
time Tx — cr^^^, r ~ 0.6332, between two different time evolution regimes. 
This crossover occurs when the noise fluctuations begin suppressing the fine 
structure of the attractor as displayed by the superstable orbit with xq = 
described previously. For r < the fluctuations are smaller than the 
distances between the neighboring subsequence positions of the xq = orbit 
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Figure 3: Two-time correlation function c{t2 — ti) for an ensemble of tra- 
jectories with Xo = for different values of noise amplitude a. See text for 
details. 

at yUc(O), and the iterate position with a > Q falls within a small band around 
the o" = position for that r. The bands for successive times do not overlap. 
Time evolution follows a subsequence pattern close to that in the noiseless 
case. When r ~ r^; the width of the noise-generated band reached at time 
Tx — 2^^"'^ matches the distance between adjacent positions, and this implies 
a cutoff in the progress along the position subsequences. At longer times 
T > Tx the orbits no longer trace the precise period-doubling structure of the 
attractor. The iterates now follow increasingly chaotic trajectories as bands 
merge with time. This is the dynamical image - observed along the time 
evolution for the orbits of a single state ^xd'^) - of the static bifurcation gap 
initially described in terms of the variation of the control parameter ji [15]. 

Two-step relaxation. Amongst the main dynamical properties displayed 
by supercooled liquids on approach to glass formation is the growth of a 
plateau, and for that reason a two-step process of relaxation, in the time 
evolution of two-time correlations [16]. This consists of a primary power- 
law decay in time difference At (so-called j3 relaxation) that leads into the 
plateau, the duration tx — t^ — I which diverges also as a power law of the 
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difference T ~Tg us the temperature T decreases to a glass temperature Tg. 
After tx tliere is a secondary power law decay (so-called a relaxation) away 
from tfie plateau [16]. In Fig. 3 we show [17] the behavior of the correlation 
function 



c(t2 - h) = 



(8) 



XhXt2 

for different values of noise amplitude. Above, (...) represents an average over 

an ensemble of trajectories starting at = and Xu = \/ — (xti)'^- The 
development of the two power-law relaxation regimes and their intermediate 
plateau can be clearly appreciated. See Ref. [10] for the interpretation of the 
map analogs of the a and f3 relaxation processes. 

Aging scaling. A second important (nonequilibrium) dynamical property 
of glasses is the loss of time translation invariance observed for T < T^, a 
characteristic known as aging. The drop time of relaxation functions and 
correlations display a scaling dependence on the ratio t/t^ where is a 
waiting time. In Fig. 4a we show [17] the correlation function 
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Figure 4: a) Two-time correlation function 
c{t + tw, tw) for different values of a. b) The same data in terms of the 
rescaled variable t/tw See text for details.] 
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C{t + ty,, ty,) = (l/N) X(t+t^)jXtj 



(9) 
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for different values of a, and in Fig. 4b the same data where the rescaled 
variable t/tw = 2" — 1, = 2A; + 1, /c = 0,1,..., has been used. The 
characteristic aging scahng behavior is patent. See Ref. [10] for an analytical 
description of the built-in aging properties of the trajectories at /Xc(c)- 

Adam-Gibbs relation. A third notable property is that the experimen- 
tally observed relaxation behavior of supercooled liquids is well described, 
via standard heat capacity assumptions [16], by the so-called Adam-Gibbs 
equation, tx — Aex.-p{B /TSc), where tx is the relaxation time at T, and the 
configurational entropy Sc is related to the number of minima of the fluid's 
potential energy surface [16]. See Ref. [10] for the derivation of the analog 
expression for the nonlinear map. Instead of the exponential Adam-Gibbs 
equation, this expression turned out to have the power law form 

tx = {s/S^f-^'^/^. (10) 

Since (1 — r)/r 0.5792 then ij, — > oo and -S'c — > as cr — > 0. 

Subdiffusion and arrest. A fourth distinctive property of supercooled 
liquids on approach to vitrification is the progression from normal diffusivity 
to subdiffusive behavior and finally to a halt in the growth of the molecular 
mean square displacement. To investigate this aspect of vitrification in the 




Figure 5: a) Repeated-cell map and trajectory, b) Mean square displacement 
{xf) for trajectories with Xq = for several values of noise amplitude a. See 
text for details. 

map at iJ^c{<^), we constructed [17] a periodic map with repeated cells of the 
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form xt+i = F{xt), F{1 + x) =1 + F{x), I = 
where 



- 1,0,1,..., F{-x) = F{x), 



, . _ J -ll-^.x^l+xa, -1 <a;<0, , . 

\1- i^,x^\+x(T, 0<x<l. 

Fig. 5a shows this map together with a portion of one of its trajectories, 
while Fig. 5b shows the mean square displacement (xf) as obtained from an 
ensemble of trajectories with xq = for several values of noise amplitude. 
The progression from normal diffusion to subdiffusion and to final arrest can 
be plainly observed as cr — > [17]. 



4 Summary 

We reviewed recent understanding on the dynamics at the onset of chaos in 
the logistic map. We exhibited links between previous developments, such as 
Feigenbaum's a function, Mori's g-phase transitions and the noise-induced 
bifurcation gap, with more recent advances, such as g-exponential sensitivity 
to initial conditions [3, 5], gr-generalized Pesin identity [4, 7] and dynamics 
of glass formation [10]. 

An important finding is that the dynamics is constituted by an infinite 
family of Mori's g-phase transitions, each associated to orbits that have com- 
mon starting and finishing positions located at specific regions of the attrac- 
tor. Thus, the special values for the Tsallis entropic index q in and Sq are 
equal to the special values of the variable q at which the g-phase transitions 
take place. 

As described, the dynamics of noise-perturbed logistic maps at the chaos 
threshold presents the characteristic features of glassy dynamics observed 
in supercooled liquids. The limit of vanishing noise amplitude o" — > (the 
counterpart of the limit T — — ^ in the supercooled liquid) leads to loss 
of ergodicity. This nonergodic state with Ai = corresponds to the limiting 
state, cr — > 0, tj; — >^ oo, of a family of small a states with glassy properties, 
which are expressed for t < tx via the g-exponentials of the Tsallis formalism. 
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